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We study the superconducting phases of the two-dimensional honeycomb lattice of graphene. We 
find two spin singlet pairing states, s-wave and an exotic p + ip that is possible because of the special 
structure of the honeycomb lattice. At half filling, the p + ip phase is gapless and superconductivity 
is a hidden order. We discuss the possibility of a superconducting state in metal coated graphene. 
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Graphene is a two-dimensional (2D) electronic sys- 
tem on a honeycomb lattice whose electronic excitations 
can be described in terms of linearly dispersing Dirac 
fermions [l| ■ Because of its unusual properties [2] , such 
as anomalous integer quantum Hall effect 0, 0] and uni- 
versal conductivity [jj], graphene has attracted a lot of 
attention in the condensed matter community. One of 
the interesting properties of graphene is that its chemi- 
cal potential can be tuned through a electric field effect, 
and hence it is possible to change the type of carriers 
(electrons or holes) , opening the doors for a carbon based 
electronics. Superconductivity has been induced in short 
graphene samples through proximity effect with super- 
conducting contacts 0. This indicates that Cooper pairs 
can propagate coherently in graphene. From the theoret- 
ical side, anomalous Andreev reflection j(| and transport 
Q have been predicted in graphene junctions with su- 
perconductors. It is known that the electronic proper- 
ties of graphene are modified by changing the number of 
graphene planes [g|. In particular, bilayer graphene has 
been demonstrated to be a tunable gap semiconductor 
0| . These results raise the question of whether it would 
be possible to modify graphene, either structurally or 
chemically, so that it would become a magnet [10( or 
even an intrinsic superconductor. By exploring the num- 
ber of graphene layers and the chemical composition, it 
maybe possible to tailor its electronic properties. 

In this letter, we derive a mean-field phase diagram 
for spin singlet superconductivity in graphene. We show 
that besides the usual s-wave pairing, an unexpected spin 
singlet state with p + ip orbital symmetry is possible 
because of the structure of the honeycomb lattice. In 
fact, we show that the p + ip state is preferred to the 
s-wave state if the on-site electron-electron interactions 
are repulsive, p + ?p-pairing states are rather interest- 
ing [ll| because, in the presence of a magnetic field, the 
superconducting vortices are described in terms of Ma- 
jorana fermions with non-abelian statistics that can be 
used for topological quantum computing Although 
it is very hard to predict pairing mechanisms from mi- 
croscopic models, we examine possible phonon and plas- 
mon mediated superconductivity in chemically modified, 
metal coated, graphene, as shown in Fig. [TJ Our results 
indicate that a plasmon mediated superconductivity is 




Figure 1: (color online) Graphene coated with metal. Small 
black (white) circles belong to the A (B) sublattices; Large 
(red) circles represent the metal atoms. 



possible in this system. 

Notice that on symmetry grounds an electronic state 
on a honeycomb lattice has three different components 
associated with the continuous SU(2) spin symmetry and 
the discrete lattice symmetry (the honeycomb lattice can 
be described as a triangular lattice with a basis with two 
atoms, A and B, as in Fig. [J). At low energies and long 
wavelengths, close to the Dirac point, the system has 
effective SO (3) Z2 spatial symmetry. Hence, the su- 
perconducting state can have a Cooper pair wavefunction 
of the form ^pair = ips ® ipz, ® i>A-B , where tps is the 
spin, tpL is the orbital, and ipA-B is the sub-lattice com- 
ponent. Pauli's principle requires the wavefunction to be 
anti-symmetric for the exchange of particles. For a spin 
singlet state, 5 = 0, one can have either L = even and 
A — B symmetric (L = being the s-wave) , or L = odd 
and A — B anti-symmetric (L = 1 being p-wave) . 

The free electron Hamiltonian can be written as: 



H, 
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where t « 2.8 eV is the hopping energy between nearest 
neighbor C atoms, a, lS (at s ) is the on-site annihilation 
(creation) operator for electrons in the sublattice A with 
spin s =T)j, and 6j iS (b\ s ) for sublattice B, n 9l , is the 
on-site particle density operator, and fi is the graphene 
chemical potential (we use units such that h = 1 = fes)- 
Diagonalization of (U) leads to an spectrum given by: 
£k = — t|7k|) where k is the 2D momentum, and 7k = 
Ea-e** = a(x/2 + V3/2y), S 2 = a(x/2 - y/3/2{}), 
and S3 — —ax, where a w 1.42 A is the C-C distance). 




Figure 2: (color online) p + ip order parameter in momentum 
space: (a) close to the Dirac point and (b) away from it. Dark 
(black) line is the real part, gray (red) is the imaginary part, 
and dashed is the amplitude of the order parameter. 

At the corners of the hexagonal Brillouin zone (at Qo = 
[0, ±47r/(3v3a)]), the band has the shape of a Dirac cone: 
£Q +k = ±t>o|k|, where vo — 3at/2 s» 6 eV A is the 
Fermi-Dirac velocity. In neutral graphene, the chemical 
potential crosses exactly through the Dirac point (// = 0). 
The electron-electron interactions are described by: 

Hp = y J2 [ a is a i* a i- s a i-s J rb\ s b is b\_ s b l - S 

is 

+5! H a \s a is h ] S ' b 0s' , (2) 
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where go and g± are on-site and nearest neighbor electron- 
electron interaction energies, respectively. It is easy to 
see that the superconducting order parameters for spin 
singlet are: (1) s-wave: Ao = (a^a^) = (hibii); (2) 
p-wave: Ai t ij = (a^bj^ — a^bji) . We assume Ai^j 
A i for all nearest neighbors and zero, otherwise. In the 
momentum space one has: Ak = JAj e~ lk '( ri ~ r j) = 
Ai7^ . Close to the Dirac points Qo (Fig.EJa)) the order 
parameter can be written as, Acj 0+ k = (3a/2)Ai(k y + 
ik x ), that is, it has p + ip symmetry. At high energies, 
away from the Dirac point, the discrete symmetry of the 
lattice is recovered and the pairing state is modified as in 
Fig. [21(b). For simplicity, however, we refer to this phase 
as p + ip, which is the symmetry close the Dirac cone. 
Decoupling the interaction terms in J2j) gives, 

Hp = E Q + g A (4i a li + b h b ll) + hx - 

i 

+9i E A Mi ( a h b ll - a l b h) + hx - > ( 3 ) 

and the total Hamiltonian can be diagonalized via a Bo- 
goliubov transformation: H cS = Y,k, a , s UJ ^a S n^ as + E , 
where E = — ,goA§ — SgiAf (A and Ai are real num- 
bers), and a is the quasi-particle number operator. 
The spectrum is: a;k, Q , s = aw^, with a, s = ±1 and, 

Wk,s = \J (i|7k| + sfif + (50A0 + . S5l A 1 | 7k |) 2 . (4) 

For Ao ^ and Ai = 0, eq. (j4]) describes an s-wave 
state with a gap given by: Eg ^ = 2|<7qAo|. In the case, 



Ao = and Ai 7^ 0, the system has an isotropic gap 
e\^ — 2\ngi Ai |/ \/t 2 + g'lA'i , which scales linearly with 
fi. In the neutral limit (p. — 0) for A = 0, the dispersion 
(jl]) is gapless, with uj^ s = i|7k|, where t is the effective 
hopping energy i — ty/l + g\A\/t 2 , which renormalizes 
the Fermi-Dirac velocity. Notice that in this case, the 
superconducting state does not lead to a gap in the spec- 
trum but to a renormalization of the velocity. This state 
of affairs we call hidden order. In the case for Ao 7^ 
and Ai 7^ 0, particle-hole symmetry is broken and the 
gap is given by E^°' 1} = 2\tg A - 0i/xAi 

The values of Ao and Ai are calculated by minimizing 
the free energy: F = £k, Q , s In (l + e"^" 3 ) + Eq, 
where (3 — l/T. The coupled self-consistent equations 
for the order parameters are: 

A = -E(goA + s|7k|5iAi)tanh(/3wk s /2)/(2w ks ), 

k,s 

A i =-£ bk| (5iAi|7k| + s5oAo)tanh(/?w k;! /2) /(6w ks ) . 

k..s 

For /i^0, one finds three distinct phases (see Fig. [3] (a)): 
(i) a s-wave phase for attractive on-site (go < 0) and re- 
pulsive nearest neighbor (<7i > 0) interactions; («) ap+ip 
phase for repulsive on-site (go > 0) and attractive nearest 
neighbor interactions (51 < 0); and (Hi) a co-existence 
phase for fully attractive interactions (go, gi < 0). The 
superconducting transitions from normal to s-wave, and 
normal to p + ip are of second order. The transitions 
involving the mixed phase for go < and <?i — > or 
gi < and go — > are of first-order even at T = 0. 
At the critical temperature, T c , the phase transitions of 
phases (i) and (it) to the normal state are second order, 
while the transition between the mixed phase and the 
normal phase is abrupt. In the weak coupling limit of 
phase (i), i.e. for \go\ <§C —go = ^Vq/K, where A is a 
high energy cut-off, the critical temperature is given by 
Q: T c « 2/i(7/7r) exp{-A( 5 g/ ff0 - l)^ 1 - 1}, where 
In 7 ~ 0.577 is the Euler constant. In reality, however, 
because the system is 2D there can be no true super- 
conducting long-range order (Mermin- Wagner theorem) 
but there will be a Kosterlitz-Thouless (KT) transition 
below a certain temperature Tkt < T c - Hence, T c only 
establishes the temperature below which the amplitude 
of the order parameter becomes finite while its phase still 
fluctuates. Phase coherence only occurs at low temper- 
atures and depends on the phase stiffness of the system. 
Although the transition is of the KT type, one expects a 
precipitous drop of the resistivity of the material below 
Tkt indicating the entrance of the electrons into a state 
of quasi-long-range superconducting order. 

For /j, = superconductivity requires a minimum cou- 
pling to occur (the problem becomes quantum critical). 
The quantum critical lines are given by go = 3§ = 
— itvq/A and g\ = g\ = — 47ri;Q/(a 2 A 3 ), as shown in Fig. 
[3] (b). We identify the phases: (iv) s-wave for go < gfi 
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Figure 3: Mean-field phase diagram: (a) /i / 0; (b) and fJ, = 0. 
Dashed lines are first order transitions, continuous lines are 
second order. 



and gi > f(g ) = 3o5i/[3o ~ (3/2) (go - .g§) 2 ]; (») gapless 
phase, for gi < h(g ) = A[v 2 F / (a 2 A 2 )]g < g\\ and (vi) 
mixed phase, for h(go) < gi < /(go), where the symme- 
try is mixed between s and p + ip pairing. 

We remark that the physical realization of supercon- 
ductivity in neutral graphene is difficult because of the 
vanishing density of states and the absence of electron- 
electron screening. Therefore, in order to superconduc- 
tivity to develop easily one has to substantially shift the 
graphene chemical potential away from the Dirac point. 
This can be achieved by chemically doping graphene with 
a metal coating: when an alkaline metal is placed on 
top of a graphene crystal, the s electrons migrate to the 
7r-band to compensate the strong difference in electro- 
negativities, raising up the chemical potential from the 
Dirac points and lowering the energy of the metallic 
bands in order to establish electrostatic equilibrium. 
Since neither carbon nor alkaline metals alone supercon- 
duct in ordinary conditions, we identify two main possible 
mechanisms for superconductivity on coated graphene: 
1) electron-phonon coupling of the ions with the electrons 
in the modified metallic bands, and 2) electron-plasmon 
coupling of the graphene electrons with the acoustic plas- 
mons of the metal. 

The electron-phonon mechanism tends to favor super- 
conductivity at high electronic densities, and has been 
used to explain the bulk superconductivity of graphite 
intercalated CaC 6 [3, [H, 0]. The strength of the 
electron-phonon coupling can be extracted from ab initio 
calculations and from experimental data, which are not 
currently available for coated graphene. Nevertheless, 
this mechanism is conventional and we will not discuss it 
here. 

In the electron-plasmon mechanism, the attractive 
electron-electron interaction is mediated by a screened 
acoustic plasmon of the metal, which is favorable to su- 
perconductivity at low electronic densities, where the 
plasmon is weakly damped by the graphene particle-hole 
continuum. In the remaining of the letter we investigate 
under what conditions this mechanism could be effec- 
tive for graphene superconductivity. Since electrons and 



acoustic plasmons have comparable energy scales, a reli- 
able calculation of the critical temperature requires the 
study of retardation effects that we will cover in future 
publications. 

Plasmon mediated superconductivity has been widely 
studied in the past 17, 3, 3, Hi 21 | and we concentrate 
on the particular aspects of the graphene problem. In 
coated graphene, the Coulomb interaction between the 
layers induces an effective electron-electron interaction 
for electrons in the graphene layer that can be calculated 
with the use of the random phase approximation (RPA) 
The RPA expansion in terms of the zeroth order 
polarization functions of the metal, n^, and of graphene, 
11°, results in an effective retarded interaction of the form: 



H ef =*52^2 V ef(<l>u)h g (ci,U})n g (-Cl,-u) , 



(5) 



where, 

Ke/(q,w)=V , g /e(q,w) [l- (Vb >q -Ki, q X(q, u)] , (6) 
is the effective electron-electron interaction, and 

e (k, w ) = i-y 0l9 [nS(k,«) + nS l (k,w)] 

+(K ) 2 9 -^ 2 9 )C(k,c)n«(k,^) (7) 

is the total dielectric function of the system, with V^q = 
2ire 2 e~ qd / (eoq) as the Fourier transform of the Coulomb 
interaction between electrons in two layers separated by 
a distance d. In Eq. ([7]), the separation between the 
metal and graphene layers induces a cross polarization 
term between the two layers that vanishes when d = 0. 
In the opposite limit, kd ^> 1, the two layers decouple. 

The electronic susceptibility of metals is commonly de- 
scribed by the Lindhard polarization function. When 
uj > vpq, where vp is the Fermi velocity of the metal, 
the q — > limit of the 2D Lindhard function gives ji^l 
Vb^II^fow) = MlJuj 2 , where Q m (q) = e v / 2E F q/e is 
the plasmon of the 2D electron gas, and Ep is the Fermi 
energy of the metallic band. The polarization function 
of graphene at small momentum q is dominated by intra- 
band excitations connecting states in the same branch of 
the Dirac cone. At lowest order, it has the same momen- 
tum and frequency dependence of the polarization func- 
tion of an infinite stack of graphite layers in the absence 
of interplane hopping 



n°(g,u;) 



-(2^/ttv 2 ) 



LU / \J~tJ 2 



vlq 2 



(8) 



In the low frequency limit, u> <C voq, the polarization 
function of graphene is approximated by its static part, 
H g (q,0), and the total dielectric function ([Tj) becomes 
e(q,uj) = e g (q,0) [l - Vl 2 p (q)/u} 2 ] , where, 



n 2 p (q)^n 2 m (q)/e g (q,0) [l- (l-e" 2 ^) V 0t9 U°(q, 0)] 



(9) 
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is the screened plasmon mode of the metal (i.e., 
e(q,Q p (q)) = 0). When qd <C 1 and \uj\ <§C v q, this 
plasmon has a linear dispersion, Sl p (q) « y / E F /(2fi)v q, 
where = [l + (8e 2 d/i/e ^)] . For > 1, the 
plasmon is not screened by the graphene layer, and one 
recovers the plasmon dispersion for the 2D electron gas: 



length limit, gii < 1, the condition for the existence of 
this acoustic mode is Ep < 2//eo/(eo + Sctgd/i/vo) < 
mvl/2 w 2.4 eV. For n ~ 2eV, e - 5 and d ~ 3 A [H, 
the left hand side of the inequality becomes E F < 1 eV, 
or equivalently that the electronic concentration in the 
metal layer has to be smaller than a c « 4 x 10 14 elec- 



Clp — ► fi TO oc In the region v F q < \u>\ -C vo<7 ; eq, ([6]) trons cm 2 (or 0.12 electrons per C). If x is the number 



can be approximated as: 

V ef (q,u,)ttV Je g (q,0)[n 2 p /(u; 2 



i] , (io) 



In the attractive region 1 of Fig. [4j the electrons of 
graphene screen the charge fluctuations, restoring the 
longitudinal response of the collective modes in the nor- 
mal phase, that is required to preserve the local gauge 
invariance of the superconductor in the London limit 
[2~3 |. The electrons of the metal by their turn are slow 
and tend to anti-screen the graphene electrons [25j, pro- 
ducing an average repulsive interaction in the metal- 
lic band. In the high frequency limit |w| 3> voq, the 
effective interaction J6]) has a second region of attrac- 
tion (region 2 of Fig. |4|). For qd <C 1, the effec- 
tive interaction in this region can be approximated by 
V e f(q,\u)\ > v q) = V , q [nj D /(uj 2 - flj D ) + l] , where 
^2_d(?) = ey/2(Ep + (J-)q/cQ is the plasmon dispersion of 
the 2D electron gas. Although the interaction is attrac- 
tive for \u>\ < D,2D(q), this attractive region is present in 
the ordinary 2D electron gas and it does not lead itself 
to superconductivity, although it favors superconductiv- 
ity by reducing the Coulomb repulsion. 

A necessary condition for appearance of plasmon in- 
duced superconductivity is that the screened acoustic 
plasmon is not overdamped by the particle-hole contin- 
uum in the interval Vp q < fl p < v$q. In the long wave- 
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Figure 4: a) Low energy excitations in the metal-graphene 
system. Dark and grey regions are the particle-hole con- 
tinuum of the metal and of graphene, respectively, b) Real 
(solid) and imaginary (dashed) parts of the effective interac- 
tion ((6]), in units of 27re 2 uo/(eoA t ) ~ 45 eV A 2 vs. momentum 
normalized by h/vq, for u> = O.l/i, /i = 2 eV, Ef = 0.4 eV. 



of remaining electrons per metallic atom M, for a system 
with chemical composition [M„C m ]2D, one of the con- 
ditions is that x < 0.12(m/n) . For K coated graphene 
KCs]2D, which is known to form a stable metallic lattice 
27|, the condition is x < 0.96 electrons per K. Values 
of x ~ 0.54—0.83 were obtained by ab initio calculations 
for K adsorbed in graphite [27j. 

In conclusion, we have derived the mean-field phase 
diagram for superconductivity in a honeycomb lattice, 
where a novel singlet p+ip phase appears. We have exam- 
ined the mechanism of plasmon mediated superconduc- 
tivity in graphene with the aim of proposing new carbon 
low-dimensional systems where superconductivity can be 
observed. 
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